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ABSTRACT 

Using the properties of SDSS DR7 QSOs catalog from Shen et al., the Baldwin effect, 
its slope evolution, the underlying drive for a large sample of 35019 QSOs with reliable 
spectral analysis are investigated. We find that the Baldwin effect exists in this large 
QSOs sample, which is almost the same in 11 different rcdshift bins, up to z ~ 5. 
The slope is -0.238 by the BCES (C iv EW depends on the continuum), -0.787 by 
the BCES bisector. For 11 rcdshift-bins, there is an increasing of the Baldwin effect 
slope from z ~ 1.5 to z ~ 2.0. From z ~ 2.0 to z ~ 5.0, the slope change is not clear 
considering their uncertainties or larger redshift bins. There is a strong correlation 
between the rest-frame C iv EW and C iv-based Mbh while the relation between the 
C iv EW and Mg n-based Mbh is very weak. With the correction of C iv-based Mbh 
from the C IV blueshift relative to Mg II, we suggest that this strong correlation is 
due to the bias of the C iv-based Mbh , with respect to that from the Mg II line. 
Considering the Mg n-based Mbh , a medium strong correlation is found between the 
C IV EW and the Eddington ratio, which implies that the Eddington ratio seems to 
be a better underlying physical parameter than the central black hole mass. 

Key words: quasars: emission lines — galaxies: active — black hole physics 



1 INTRODUCTION 

Broad emission lines are a prominent property of QSOs. 
The discovery of an anti-correlation between the equiva- 
lent width (EW) of C IVA1549 emission line and its nearby 
continuum luminosity in QSOs rest frame (the Baldwin ef- 
fect), was first made by Baldwin (1977) (e.g. see a review by 
Shields 2006). Over the past 30 years, a significant amount 
of effort has been expended to confirm this effect in C iv 
including other prominent emission lines (e.g., Lya, C III , 
Si IV , Mg II, [O ill] , Fe Kq), as well as exploring its origin 
and evolution (e.g. Dietrich et al. 2002; Zhou et al. 2005; 
Netzer et al. 2006; Kong et al. 2006; Xu et al. 2008; Wu et 
al., 2009; Richards et al. 2011). 

Although it is believed that the Baldwin effect exists for 
many UV/optical emission lines, its origin is still a problem 
to debate. Several interpretations about its origin have been 
proposed (e.g. Netzer et al. 1992; Dietrich et al. 2002; Shang 
et al. 2003; Baskin & Laor, 2004; Xu et al. 2008; Wu et al., 
2009). One promising interpretation is the softening of the 
spectral energy distribution (SED) for increasing luminos- 
ity, which lowers the ion populations having high ionization 
potentials (e.g., Netzer et al. 1992; Dietrich et al. 2002). 

With the progress in the virial mass of a supermassive 
black hole (SMBH, Mbh ), it provides the possibility to ex- 



plore the underlying physical parameters for the origin of C 
iv Baldwin effect, such as the Eddington ratio (i.e. the ra- 
tio of the bolometric luminosity to the Eddington luminos- 
ity, -^Boi/i/Edd ), the Mbh , and the luminosity dependence 
of metallicity (e.g. Dietrich et al. 2002; Shang et al. 2003; 
Warner et al., 2004; Baskin & Laor, 2004, 2005; Xu et al. 
2008). The single-epoch virial Mbh can be calculated from 
the broad line width (e.g., H/3 , H« , Mg II, C iv) and empir- 
ical luminosity- size relation (e.g. Kaspi et al. 2000; Bian & 
Zhao 2004; Vestergaard & Peterson 2006; Shen et al. 2011). 
The LboI can be calculated from the monochromatic lumi- 
nosity by the correction factor estimated from the composite 
SED (e.g., Richard et al., 2006; Shen et al. 2011). 

For PG QSOs sample, Baskin & Laor (2004) found a 
strong correlation between the C iv EW and LBoi/^Edd , 
and suggested that the LBoi/^Edd is the primary physical 
parameters driving the Baldwin effect. However, Shang et 
al. (2003) used the method of Spectral Principal Component 
(SPC) to find there is no correlation between the Baldwin 
effect and the -i/Boi/^Edd or the Mbh as underlying physi- 
cal parameters. Nikolajuk & Walter (2012) found that weak 
line QSOs didn't follow the relation between C iv EW and 
iBoi/^Edd by Baskin & Laor (2004). However, these weak- 
line QSOs have the same x-ray-optical spectral index (cw), 
with respect to other normal QSOs. It was suggested that 
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the weak-line QSOs are caused by the possible low covering 
factor (Nikolajuk & Walter 2012). 

In the past, most of the work was made by using rela- 
tively small samples. The Sloan Digital Sky Survey (SDSS) 
provides possibility for us to investigate the Baldwin effect 
in a larger QSOs sample (e.g. Xu et al., 2008; Richard et al., 
2011). With SDSS DR5, Xu et al. (2008) used the result from 
the pipeline of SDSS spectral fit, which measured the line 
feature by a single Gaussian fit for the continuum subtracted 
spectrum. Xu et al. (2008) found that, up to z ~ 5, the slope 
of the Baldwin effect seems to have no effect of cosmologi- 
cal evolution. They also found the C IV EW has a stronger 
correlation with C iv-based Mbh than the I/Boi/^Edd , sug- 
gesting that the Mbh is probably the primary drive for the 
Baldwin effect. 

Recently, Shen et al. (2011) made a careful spectral 
analysis and gave a compilation of properties for 105783 
QSOs from SDSS DR7 QSOs catalog. Here we use their data 
to reinvestigate the Baldwin effect for the C IV line and its 
origin. Our adopted sample is briefly described in §2, the 
results and the analysis are given in §3, and the conclusions 
are presented in §4. 



2 SAMPLE 

SDSS DR7 covers an imaging area of about 11663 square de- 
grees and a spectroscopic area of about 9380 square degrees 
(Abazajian et al. 2009). Schneider et al. (2010) presented a 
quasars catalog from SDSS DR7, consisting of 105783 QSOs 
with Mi < —22, increasing by 30% with respect to that in 
DR5. The wavelength coverage of SDSS spectrum is 3800A 
to 9200A. In order to investigate the C IV line, we select 
51523 QSOs with z > 1.5 in SDSS DR7 QSOs. 

Shen et al. (2011) carefully considered the UV/optical 
Fe II contribution, host contribution in low-z QSOs, and 
multiple-Gaussian fit for four lines of Ha , H/3 , Mg II, C IV 
for the single-epoch virial SMBH mass calculation. In their 
SDSS spectral analysis, they remove the effects of Galactic 
extinction in the SDSS spectra, and shift the spectra to rest 
frame. They compiled these continuum and line properties, 
as well as radio, infrared , X-ray properties (Schneider et al. 
2010). 

Because adopting the properties of the rest-frame C IV 
EW, Mbh and I/Boi/^Edd in Shen et al. (2011), we briefly 
described their spectral analysis. For the weakness of Fe II 
and not high enough signal-to-noise (SN) of C IV line, Shen 
et al. (2011) only reported their C IV measurements with- 
out the Fe II template fits, and emphasized that their C IV 
EWs may be averagely overestimated by 0.05 dex. They fit 
for the C IV line over the [1500A, 1600A] wavelength range, 
with multiple Gaussians. They didn't subtract the C IV nar- 
row component (e.g., Sulentic et al. 2007). Considering the 
possible narrow or broad absorption features of the C IV line 
in high-z QSOs, they masked out 3<r outliers below the 20 
pixel boxcar-smoothed spectrum during the fits for the nar- 
row absorption, and performed a second fit excluding pixels 
below 3cr of the first model fit for the broad absorption. 
They gave the BAL (broad absorption line) flag for SDSS 
DR7 QSOs. QSOs flagged as BAL QSOs are not used in our 
investigation of the Baldwin effect. 

For QSOs with 1.5 < z < 2.25, they also fit the Mg 
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Figure 1. The distributions of the redshift z, SN(C iv), and 
reduced x 2 - 



II emission line, as well as the SMBH mass from Mg II 
line. They fit for the Mg II line over the [2700A, 2900A] 
wavelength range, with a single Gaussian (with FWHM < 
1200km s" 1 ) for the narrow component, and for the broad 
component with a single Gaussian or multiple Gaussians. 
They provided a new mass calibration from Mg II lines by 
their multiple- Gaussian fits with narrow line subtraction, 
which is adopted as the fiducial virial SMBH mass for QSOs 
with 0.9 < z < 1.9. For z > 1.9, they adopted the formulae 
of Vestergaard & Peterson (2006) to calculate the C IV virial 
mass. 

We use following criteria to select our sample: z > 1.5, 
FWHM(C iv) > 1200 km s" 1 , reduced x 2 < 2.0, EW(C iv) 
> 2.0, SN(C iv) > 5, excluding BAL QSOs. The number of 
our QSOs sample is 35019 (~ 68% of 51523 QSOs). Their 
distributions of the redshift, SN, and reduced x 2 are showed 
in Figure 1. For the SMBH mass, we adopted their fiducial 
viral masses from Mg II for z < 1.9 and masses from C IV for 
z > 1.9. The numbers are 14652, 19476, respectively (Shen 
et al. 2011). In Xu et al. (2008), there are 26623 QSOs used 
to investigate the C iv Baldwin effect, and 13960 QSOs used 
to investigate the relation of C iv EW with the accretion 
parameters (M B h , L Bo i/L Edd )■ 

3 RESULTS AND DISCUSSION 

3.1 The relation between the rest- frame C iv EW 

and I/i35o 

The following formula is used to express the relation between 
the C iv EW and the continuum luminosity at 1350 A (L1350 
). 

logEW(C IV) =Q+/31ogLi350 (1) 

Figure 2 gives the relation between the C iv EW and L1350 
for 35019 QSOs from SDSS DR7. For the total 35019 QSOs, 
the spearman's rank correlation coefficient is R — —0.33 and 
the null hypothesis is less then 10~ 4 . 

The regression analyzes are performed using the bivari- 
ate correlated errors and intrinsic scatter (BCES) algorithm 
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Figure 2. log-log correlation between the C IV EW and the con- 
tinuum luminosity at 1350A (in erg s — 1 ) for the DR7 sample of 
35019 QSOs with 1.5 < z < 5.0. In the left bottom corner, the 
numbers for colorful contour levels are shown. The solid line is the 
linear fit from the BCES bisector and the bootstrap simulation. 
The dash line is BCES with EW(C iv) as the dependent variable 
and Z/1350 as the independent variable. In the right top corner, 
the redshift coverage, the number of the sample, the spearman's 
rank correlation coefficient, and the slopes with BCES(y|X) and 
BCES bisector are listed. Black circles arc data from a compiled 
sample of Wu et al. (2009). 



(Akritas & Bershady 1996) . There exists some difference be- 
tween the BCES(Y \X) [i.e., Y = f(X)] and BCES(Y| Y) [i.e., 
X = g(Y)] regressions. The bisector bisects these two regres- 
sions. The BCES regression is quite robust; bootstrapping 
simulations reproduce the theoretically expected results well 
(see, e.g., Akritas & Bershady 1996). For the SDSS sam- 
ple, with the BCES(Y|Y) (X is Li 3S o ), the slope is j3 = 
-0.238+0.040, and the intercept is a = 12.6T0.19. However, 
with the BCES bisector, the slope is /3 = -0.787 ± 0.026, 
and the intercept is a = 37.9 ± 1.20 (see Table 1). 

In Figure 2, we also plot their compiled sample from Wu 
et al. (2009). Their compiled sample consists of 272 AGNs: 
189 from SDSS, 47 from HST, and 36 from IUE. They used 
two Gaussian profiles to fit the C IV emission lines after 
subtracting the power-law continuum. The continuum lumi- 
nosity at the rest-frame 2500A in Wu et al. (2009) is con- 
verted to the luminosity at the rest-frame 1350A by assum- 
ing /„ oc i>~ ' 44 (i.e, larger by ~ 0.12dex, Vanden Berk et 
al. 2001). In Figure 2, the data of Wu et al. (2009) are con- 
sistent with our SDSS DR7 sample, as well as the dash line 
from BCES(Y|Y). For their compiled sample of Wu et al. 
(2009), there are 163 common QSOs in our SDSS DR7 sam- 
ple. For these 163 common QSOs, the C iv EWs from Wu et 
al. (2009) are consistent with that from Shen et al. (2009) 
very well, their C iv EW ratio has a normal distribution 
with mean 0.01 and dispersion 0.12. 

For 81 PG QSOs from Baskin & Laor (2004), 454 Large 
Bright Quasar Survey objects by Forster et al. (2001) and 
125 pre-COSTAR AGN by Kuraszkiewicz et al. (2002), we 
also find that these QSOs/ AGNs are located in the re- 
gion of SDSS QSOs or consistent with the dash line from 
BCES(Y|X) if they are scaled to our Figure 2. 

The Baldwin effect provides a potential method to infer 
the luminosity of a QSO from its C iv observation. With 



respect to classical standard candles of type la supernovae 
(SNe), QSOs can be observed to a much higher redshift. For 
our SDSS DR7 sample, we find the relation between L1350 
and rest-frame C iv EW by BCES (Y\X), BCES bisector, 
respectively: 

logLiaso = (-0.511 ± 0.042) logEW(C iv) 
+(46.977 ± 0.066) 

log £1350 = (-1.271 ± 0.043) logEW(C iv) 
+(48.179 ±0.068) 

From above formula, we can calculate the rms values of the 
residuals after subtracting the predicted luminosity from the 
measured C iv EW. We find that they are 0.286, 0.345, 
respectively. In order to use QSOs as standard candles via 
the Baldwin effect, it is suggested that we should at least 
confine the luminosity within an uncertainty of 30%, ~ 0.13 
dex (Wu et al., 2009). Therefore, it is impossible to use our 
QSOs data-set as standard candles via the Baldwin effect. 

3.2 The slope of the Baldwin effect 

In previous studies, the slope of the Baldwin effect is 
-0.14 r 0.20 (see Table 1 in Xu et al. 2008 and refer- 
ence therein). Dietrich et al. (2002) found the slope becomes 
steep for luminous QSOs , from —0.14 for their total sam- 
ple to -0.20 for QSOs with log AL a (1450 A) > 44 erg s" 1 . 
We get a slope from the BCES(Y|X), p = -0.238 ± 0.04 
(Table 1), which is steeper than some previous studies. By 
BCES bisector, the slope is /? = -0.787 ± 0.026. The slope 
of BCES bisector (—0.787) is more consistent with the slope 
of the "intrinsic" Baldwin effect, (/3 = -0.72) by Pogge & 
Peterson (1992). Considering the compiled sample from Wu 
et al. (2009), the slope from the BCES(Y|Y) is better than 
the slope of BCES bisector when wide luminosity range for 
QSOs are concerned (Figure 2). 

The redshift coverage of our sample is from 1.5 to 5.0 
(see Figure 1). In order to investigate the relation of the 
Baldwin effect with the redshift, We divide our sample into 
11 bins according to the redshift, and the redshift bins are 
1.5-1.6, 1.6-1.7, 1.7-1.8, 1.8-1.9, 1.9-2.0, 2.0-2.1, 2.1-2.25, 
2.25-2.5, 2.5-3.0, 3.0-3.5, 3.5-5.0 (Table 1). The number in 
different redshift bins is about 3000. For different redshift 
bins, we calculate the spearman's rank correlation coeffi- 
cient, the null possibility, the slopes and the intercepts with 
the BCES (Y\X) and the BCES bisector (Table 1). Except 
the highest bin (3.5 < z < 5.0), We find a medium strong 
correlation for different redshift bins (\R\ > 0.3). 

It seems that there is an increase of the slope from z ~ 
1.5 to z ~ 2 (Figure 3), which is the same to that by the 
BCES (Y\X) and BCES bisector. From z ~ 2 to z ~ 5, the 
slope evolution is not clear due to the large uncertainties. If 
we believe that the Baldwin effect is due to the softening of 
SED for increasing luminosity, the slope evolution from z ~ 
1.5 to z ~ 2 implies the evolution of QSOs SED below z ~ 2. 
It is also possibly related to the covering factor evolution. In 
low redshift bins, the dynamical range in luminosity is larger 
due to the flux-limited nature of the sample, which might 
affect the linear regression. For redshift bins below z — 2.0, 
in the same luminosity range, i.e., 10 45 ' 6 to 10 46 ' 6 erg s _1 , 
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Table 1. The C IV EW and the continuum luminosity Regression Parameters. logEW(C iv) = a + /31ogLi35o. 



z 


N 


R 


/3 1 


Q 1 


/? 2 


a 2 


-fnull 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 




SDSS DR7 sample 


1.5-5.0 


35019 


-0.33 


-0.238 ±0.040 


12.6 ±0.19 


-0.787 ± 0.026 


37.9 ± 1.20 


1.15 x 10" 


-25 


SDSS DR7 sample in different z bins 


1.5-1.6 


3699 


-0.33 


-0.386 ±0.021 


19.4 ±0.98 


-1.065 ± 0.015 


50.2 ±0.70 


3.90 x 10 


-4 


1.6-1.7 


3645 


-0.32 


-0.292 ±0.016 


15.0 ±0.75 


-0.911 ± 0.020 


43.5 ±0.93 


8.93 x 10 


-4 


1.7-1.8 


3727 


-0.35 


-0.283 ±0.013 


14.6 ±0.58 


-0.887 ±0.014 


42.4 ±0.63 


2.65 x 10" 


-5 


1.8-1.9 


4057 


-0.35 


-0.272 ±0.012 


14.1 ±0.57 


-0.845 ± 0.014 


40.5 ±0.65 


5.29 x 10 


-5 


1.9-2.0 


3343 


-0.31 


-0.233 ±0.013 


12.3 ±0.60 


-0.830 ±0.014 


39.9 ±0.63 


2.73 x 10 


-3 


2.0-2.1 


2687 


-0.32 


-0.230 ±0.015 


12.2 ±0.70 


-0.858 ± 0.015 


41.2 ±0.67 


6.10 x 10 


-3 


2.1-2.25 


3265 


-0.33 


-0.249 ±0.013 


13.1 ±0.62 


-0.850 ±0.016 


40.9 ±0.74 


1.32 x 10 


-3 


2.25-2.5 


3173 


-0.34 


-0.273 ±0.014 


14.2 ±0.64 


-0.840 ±0.015 


40.4 ±0.67 


4.99 x 10 


-4 


2.5-3.0 


2740 


-0.34 


-0.279 ±0.015 


14.5 ±0.69 


-0.845 ± 0.023 


40.7 ±1.06 


9.22 x 10 


-4 


3.0-3.5 


3213 


-0.30 


-0.238 ±0.014 


12.6 ±0.66 


-0.698 ± 0.242 


33.9 ±11.2 


7.52 x 10 


-3 


3.5-5.0 


1470 


-0.27 


-0.215 ±0.022 


11.5 ± 1.04 


-0.031 ± 0.696 


2.93 ±32.4 


1.59 x 10 


-1 






SDSS DR7 sample in different z bins (10 45 


6 < ^1350 < 10 


46.6gj.g3-l) 






1.5-1.6 


3566 


-0.30 


-0.387 ±0.023 


19.4 ± 1.05 


-1.140 ±0.015 


53.8 ±0.70 


4.74 x 10" 


-3 


1.6-1.7 


3474 


-0.28 


-0.293 ±0.018 


14.9 ±0.83 


-0.997 ± 0.018 


47.5 ±0.81 


1.43 x 10" 


-2 


1.7-1.8 


3836 


-0.31 


-0.287 ±0.016 


14.7 ±0.74 


-0.919 ± 0.012 


44.1 ±0.55 


1.41 x 10" 


-3 


1.8-1.9 


3518 


-0.31 


-0.280 ±0.016 


14.5 ±0.73 


-0.882 ± 0.052 


37.6 ±0.24 


3.52 x 10 


-3 


1.9-2.0 


3136 


-0.27 


-0.230 ±0.016 


12.2 ±0.75 


-0.828 ± 0.011 


38.9 ±0.53 


3.48 x 10 


-12 



1 :BCES(y|X) where Y is EW(C iv) as the dependent variable, and X is L1350 as the independent variable. 2 :BCES bisector result. 



we find R is ~ —0.3, and the slopes are consistent with the 
results for the total sample in different redshift bins (Table 

1)- 

With SDSS DR5, Xu et al. (2008) used the SpecLine 
table to investigate C iv Baldwin effect. The SpecLine table 
is from the pipeline of SDSS spectral fit. The C iv line fea- 
ture is measured by a single Gaussian fit for the continuum- 
subtracted spectrum. Xu et al. (2008) investigated the C iv 
Baldwin effect in the observed frame. In Xu et al. (2008), 
there was no correction of Galactic extinction for the L1350 
and the C IV EW was given in the observed frame. In Fig- 
ure 4, we showed a comparison of C iv EW in the rest frame 
for 10232 common QSOs from the current sample and Xu 
et al. (2008). We find that, the rest-frame C iv EW from 
the current sample is averagely larger by 0.23 dex than that 
from Xu et al. (2008). The difference is due to the spectral 
fitting, such as a single Gaussian fit for C iv in Xu et al. 
(2008), multiple-Gaussian fit in Shen et al. (2011). 

3.3 The relation between C iv EW and A/bh or 

iBol/^Edd 

The underlying physical origin for the Baldwin effect is still 
an open question. The fundamental process in QSOs is the 
accretion around the SMBHs. Using the data from Shen et 
al. (2011), we also investigate the relation between the rest- 
frame C iv EW and Mbh or Lsoi/iEdd (Figure 5). For the 
total sample, we find a medium strong correlation between 
the rest-frame C iv EW and M B h with R = -0.27, and 
the null hypothesis is less then 10~ 4 . However, the corre- 
lation between the C iv EW and Z/Boi/^Edd is very weak, 
R = 0.01 (Figure 5). Earlier studies on individual sources 
suggest that for the same object, an intrinsic Baldwin effect 



does exist (albeit with a different slope), therefore, luminos- 
ity does drive the Baldwin effect (for the same object, it is 
equivalent to say that I/Boi/iEdd drives the Baldwin effect). 
We find that the very weak correlation between C iv EW 
and LBoi/^Edd for the total SDSS DR7 sample is due to the 
mixed use of Mg 11-based and C iv-based Mbh estimates. 

For QSOs with z > 1.9, we use the mass from the C 
iv line, while the mass from Mg 11 line for QSOs with 1.5 < 
z < 1.9 (Shen et al. 2011). So, we divide the total our sample 
into two samples, i.e., 1.5 < z < 1.9 and z > 1.9. For QSOs 
with 1.5 < z < 1.9, We find a very weak correlation of 
R = 0.07 between the rest-frame C iv EW and mass from 
the Mg 11 line. For the physical drive of the Baldwin effect, 
L1350 and Mbh are two essentially independent quantities. 
To understand which one is a better physical drive on the 
Baldwin effect, one needs to look at them separately, i.e., 
at fixed Z/1350 or fixed Mbh ■ Using the whole sample with 
mixed Mbh and L1350 makes it difficult to isolate the two. 
For L1350 between 10 45 7 and 10 46 5 erg s _1 , we divide the 
total SDSS sample into four luminosity bins (AlogLisso = 
0.2), and find that this correlation between C ivEW and 
M B h is stronger (R ~ 0.18) than that for the total SDSS 
sample (R ~ 0.07). For QSOs with z > 1.9, we find a strong 
correlation of R = —0.48 between the C iv EW and C IV- 
based Mbh (Table 2; Figure 6). It is consistent with our 
previous result where we used the C iv FWHM from the 
standard pipeline of SDSS to calculate the single-epoch Mbh 
(Vestergaard & Peterson 2006; Xu et al. 2008). In section 
3.4, the bias in C iv-based Mbh will be discussed. 

For QSOs with 1.5 < z < 1.9, we find a modest cor- 
relation between C iv EW and L Bo i/L Edd (R = —0.28). 
Because Mg 11-based Mbh is proportional to v 2 x -Rblr, and 
Rbl,r oc L1350 (Shen et al. 2011), I/Boi/iEdd is proportional 
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Figure 3. The slope versus the redshift. The red points are for the slopes by the BCES(Y|X). The green ones are for the slopes by the 
BCES bisector. The red/green dash lines are for the slopes for the total sample, respectively by the BCES(Y|X) and bisector. 
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Figure 4. Top: A comparison of the rest-frame C IV EW for 10232 common QSOs from the current sample and Xu et al. (2008). 
dash line is 1:1. Bottom: the distribution of the EW ratio of QSOs from Xu et al. (2008) to the current sample. 
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Table 2. The C IV EW drive Regression Parameters. logEW(C iv) = a + /3 log Li 350 , logEW(C iv) = a + j3 log M B h , or logEW(C iv) = 
a + /31ogL Bo i/L E dd- 



Drive 


z 


N 


R 


/3 1 


a 1 




a 2 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


£-1350 


1.5-5.0 


35019 


-0.33 


-0.238 ± 0.040 


12.6 ±0.19 


-0.787 ± 0.026 


37.9 ± 1.20 




1.5-1.9 


15128 


-0.34 


-0.299 ± 0.076 


15.4 ±0.35 


-0.920 ± 0.007 


44.0 ± 0.34 




1.9-5.0 


19891 


-0.31 


-0.226 ±0.051 


12.0 ±0.23 


-0.755 ±0.057 


36.5 ± 2.63 


M BH 


1.5-5.0 


34128 


-0.27 


-0.125 ±0.011 


2.74 ± 0.10 


-0.731 ±0.025 


8.34 ±0.23 




1.5-1.9 


14652 


0.07 


0.665 ± 0.087 


0.985 ± 0.08 


0.959 ± 0.057 


-7.22 ±0.05 




1.9-5.0 


19476 


-0.48 


-0.178 ±0.019 


32.3 ± 0.18 


-0.521 ±0.052 


6.41 ± 0.49 


Lboi/ ^Edd 


1.5-5.0 


34128 


0.01 


0.272 ± 0.001 


1.58 ± 0.01 


0.182 ±0.091 


1.77 ±0.83 




1.5-1.9 


14652 


-0.28 


-0.268 ± 0.080 


1.41 ± 0.06 


-0.858 ± 0.080 


0.992 ± 0.07 




1.9-5.0 


19476 


0.26 


0.0003 ± 0.00009 


1.57 ±0.01 


0.090 ± 0.067 


1.69 ± 0.07 



1 :BCES(y|X) where Y is EW(C iv) as the dependent variable, and X is Z/1350 , Mbh , ^Bol/^Edd as t ne independent variable, 

respectively. 2 :BCES bisector result. 



© 0000 RAS, MNRAS 000, 000-000 



6 Bian et al. 



1.50<=z<5.00 

Num=34128 

R=-0.27 

Slope:-0.731( 0.025); - 




£ 2 



1.50<=z<5.00 
Num=34128 
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Figure 5. log-log correlation between the rest-frame C IV EW and the A/bh (left), iBol/^Edd (right)- In the right of two panels, the 
numbers for colorful contour levels are shown. The solid line is the linear fit from the BCES bisector and the bootstrap simulation. The 
dash line is BCES with EW(C iv) as the dependent variable and A/bh /-^Bol/^Edd as the independent variable. In the right top corner, 
the redshift coverage, the number of the sample, the spearman's rank correlation coefficient, the slopes and their errors with BCES(y|X) 
and BCES bisector are listed. 



t.50<-z<1.90 
Num-14652 
R= 0.07 

slope: 0.959( 0.006); 0.0Bo( 0.009) 
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Num-14652 

R=-0.28 

slope:-0.859( 0.008); -0.268( 0.008) 
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R= 0.26 

slope: 0.094( 0.066); 0.000( 0,000) 




-0.5 
'°9 <^„/kdd) 



-1 -0.5 
lo 9 (l^/kdd) 



Figure 6. log-log correlation between the rest-frame C IV EW and the A/bh (top), X/Bol/^Edd (bottom) for two redshift bins, 1.50 < 
z < 1.9, z > 1.9. In the right of four panels, the numbers for colorful contour levels are shown. The solid line is the linear fit from the 
BCES bisector and the bootstrap simulation. The dash line is BCES with EW(C iv) as the dependent variable and A/bh /^Bol/^Edd 
as the independent variable. In the right top corner of four panels, the redshift coverage, the number of the sample, the spearman's rank 
correlation coefficient, the slopes and their errors with BCES(y|X) and BCES bisector are listed. Black circles in two bottom panels are 
data from Baskin & Laor (2004). 
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to L?35o neglecting the effect of Mg II FWHM. We find that 
EW(CIV) oc (L B oi/L E dd) _0 ' 27 for QSOs with 1.5 < z < 1.9 
by BCES(Y \X). Using this correlation, we can derive the ex- 
pected Baldwin effect of EW(CTV) oc -^1350^: however, the 
slope of Baldwin effect (Table 1) is -0.24. It suggests that the 
relation between C IV EW and I/Boi/^Edd is not completely 
from the relation between C IV EW and L1350 • We use par- 
tial Kendalls t to investigate the role of L1350 in the relation 
between C iv EW and I/Boi/^Edd for the total sample. We 
find that the relation between C IV EW and LBoi/iEdd is 
weakened when L1350 is held fixed. The partial Kendalls r 
is -0.135 for QSOs with 1.5 < z < 1.9, and 0.097 for QSOs 
with 1.9 < z < 5.0, while their spearman's rank correlation 
coefficients are -0.28 and 0.26 (Table 2). It is possible that 
Lsoi/LEdd is not very important in controlling the Bald- 
win effect, or that the large uncertainty of individual Mbh 
estimates is diluting any inherent correlation. 

In Figure 6, we also show the location of low-z PG QSOs 
from Baskin & Laor (2004), where the Mbh was calculated 
from the H/3 FWHM. These PG QSOs, which are denoted by 
black circles in Figure 6, are consistent well with our sample 
of 1.5 < z < 1.9 QSOs where Mg n-based Mbh is used. 
However, these PG QSOs don't follow the relation found 
for 1.9 < z < 5.0 QSOs where C iv-based Mbh is used. 
With respect to H/3 -based I/Boi /-I/Edd , it was found that 
the relation between EW(C iv) and I/Boi/^Edd is weaker 
when LBoi/^Edd is based on the C iv FWHM (Baskin & 
Laor 2005). Considering the strong correlation between C 
iv EW and LBoi/^Edd for low-z PG QSOs found by Baskin 
& Laor (2004), our result also shows a strong bias of C IV- 
based Mbh , with respect to Mg n-based Mbh • 

3.4 The C iv-based Mbh correction from the C iv 
blueshift relative to Mg 11 

However, we should keep in mind that the underlying physi- 
cal drive analysis depends on the Mbh calculation. The bias 
of Mbh calculated from Mg n, C iv are investigated (e.g., 
Rafiee & Hall 2011; Shen et al. 2012). It is found that the C 
iv-based Mbh is biased to the possible non-virialized com- 
ponent in C iv (Shen et al., 2008; Richard et al. 2011; Shen 
et al. 2011; Shen et al. 2012). The Mg 11-based Mbh can 
be calibrated to yield consistent virial mass estimates with 
those based on the Ha /H/3 , while the C iv-based Mbh is 
poorly correlated with the H/3 -based A/bh or Ha -based 
Mbh (Shen et al. 2012). It was found that the C IV blueshift 
relative to H/3 can be used as a C iv FWHM correction 
(shen et al. 2012). For 24228 QSOs with 1.5 < 2 < 2.25, we 
find a strong correlation between C iv-based Mbh and the 
C iv blueshift relative to Mg 11 (AV = u ff,civ — v b,m s u), 
R — —0.41. We also find a strong correlation between C iv 
EW and the C iv blueshift relative to Mg 11, R = -0.38 
(bottom left panel in Figure 7; R — —0.11 when only using 
the C iv blueshift). This strong correlation between C iv 
EW and the C iv blueshift relative to Mg 11 is consistent 
with the disk+wind model (Richards et al. 2011). It is pos- 
sible that QSOs with larger C iv blueshift have softer SEDs 
and have a larger radiation line driving contribution to their 
winds to dominate over the disk. 

The C iv blueshift relative to Mg 11 would lead to the 
apparent correlation between C iv EW and C iv-based Mbh 
. We can use the C IV blueshift relative to Mg 11 to do the 



C iv-based Mbh correction. Shen et al. (2011) gave the Mg 
11-based Mbh and C iv-based Mbh for QSOs with 1.5 < 
z < 2.25 (their Table 1). Their spearman's rank correlation 
coefficient is R = 0.28 and the null hypothesis is less then 
10~ 4 (top left panel in Figure 7). We use the C IV blueshift 
relative to Mg 11 to do the correction of C iv-based Mbh • 
Using BCES regression analyzes, we find a strong correlation 
between the Mbh difference (logMsH.Mgii — logMsH.Mgii) 
and the C iv blueshift relative to Mg 11, R=-0.44 (top right 
panel in Figure 7): 

log M B H,M g n - log Mbh, civ = a + 6(u ff ,crv - «off,M g ii) (2) 

Where a = (0.174 ± 0.004), b = (-2.01 ± 0.04) x 10~ 4 by 
BCES(y \X), and a = (0.449 ± 0.007), b = (-5.78 ± 0.08) x 
10~ 4 by BCES bisector. 

Considering the correction of C iv-based Mbh from 
above formula, the spearman's rank correlation coefficient 
for the relation of Mg 11-based Mbh and the corrected C 
iv-based Mbh changes from 0.28 to 0.38, 0.36, respectively. 
Using above formula by BCES(F|A) and BCES bisector to 
do the C iv-based Mbh correction, we find that, for the re- 
lation between C iv EW and Lsoi/iEdd for 9019 QSOs with 
1.9 < z < 2.25, R changes from 0.22 to -0.03, -0.29, respec- 
tively (bottom right panel). And for the relation between 
C iv EW and the corrected C iv-based Mbh , R changes 
from -0.45 to -0.24, 0.1, respectively. They are consistent 
with the results for 14652 QSOs with 1.5 < z < 1.9 where 
Mg 11-based Mbh is used. 



3.5 The origin of the Baldwin effect 

Considering Mg 11-based Mbh (QSOs with z < 1.9 QSOs), 
we find that the correlation coefficient for the relation be- 
tween C iv EW and L1350 {R = —0.34) is larger than the 
other two relations (i.e. C iv EW versus A4bh , C iv EW 
versus LBoi/^Edd ; R = 0.07, —0.28). And the correlation co- 
efficient of C iv EW with L Bo i/L Edd (R = —0.28) is larger 
than that of C iv EW with Mbh (R = 0.07). It is consis- 
tent with the result by Baskin & Laor (2004). The stronger 
correlation between C iv EW and Mbh in Xu et al. (2008) 
is due to the usage of C iv-based Mbh • Therefore, by the 
larger sample, we find that C IV EW is primarily controlled 
by the L1350 , with respect to -I/Boi/^Edd and Mbh ■ The 
modest correlation between C iv EW and the Lboi/Lem 
(larger than that between C iv EW and Mbh ) implied that 
the I/Boi/^Edd seems to be a better underlying physical pa- 
rameter than Mbh ■ The strong correlation between the C 
iv EW and C iv blueshift relative to Mg 11 is in favor of the 
disk+wind model (e.g. Richards et al. 2011). The QSOs with 
larger C iv blueshift have softer SEDs and have a larger ra- 
diation line driving contribution to their winds to dominate 
over the disk. 



4 CONCLUSIONS 

With respect to our previous analysis for SDSS DR5 QSOs, 
we use the SDSS DR7 QSOs catalog from Shen et al. (2011) 
to reinvestigate the Baldwin effect, its slope evolution, the 
underlying drive. The main conclusions can be summarized 
as follows: (1) The Baldwin effect exists in the large sample 
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Figure 7. Top left: C iv-based Mbh versus Mg Il-based Mbh ■ Top right: The ratio of Mg Il-based Mbh to C iv-based Mbh versus 
the C IV blueshift relative to Mg II. Bottom left: C IV EW versus the C IV blueshift relative to Mg II. Bottom right: C IV EW versus 
the corrected C iv-based i/Bol/^Edd (f° r clarity, the data points are not shown). The numbers for colorful contour levels are shown in 
different panels. The solid line is the linear fit from the BCES bisector and the bootstrap simulation. The dash line is BCES with Y as 
the dependent variable and X as the independent variable. In the top corner of four panels, the redshift coverage, the number of the 
sample, the spearman's rank correlation coefficient, the slopes and their errors with BCES(Y|X) and BCES bisector are listed. Black 
circles in the bottom right panel are data from Baskin & Laor (2004). 



of 35019 QSOs with reliable spectral analysis. The correla- 
tion coefficient is R — —0.33, which is almost the same for 
QSOs in different redshifts, up to z ~ 5. (2) For the total 
sample, the slope is -0.238 by the BCES (C IV EW depends 
on L1350 ), -0.787 by the BCES bisector. (3) For 11 red- 
shift bins, there is an increasing of the Baldwin effect slope 
from z ~ 1.5 to z ~ 2.0. From z ~ 2.0 to z ~ 5.0, the 
slope change is not clear considering their uncertainties or 
the larger redshift bin. (4) By 34128 QSOs with M BH and 
Lsoi/L-Edd carefully derived from spectral decomposition by 
Shen et al. (2011), we find that there is a strong correlation 
between the rest-frame C IV EW and C iv-based Mbh for 
1.9 < z < 5.0(R = —0.48). However, the relation between 
C IV EW and Mg il-based Mbh for 1.5 < z < 1.9 is very 
weak (7? = 0.07). We think it is due to the bias of the C iv- 
based Mbh , with respect to the Mg n-based Mbh • Using 
the correction of C iv-based Mbh from the C IV blueshift 
relative to Mg II, we find the correlation of C IV EW with 
corrected C iv-based Mbh becomes weak (from R = —0.45 
to R = 0.1). Considering the Mg n-based Mbh , a medium 
strong correlation is found between C IV EW and LBoi/^Edd 



(R — —0.28), which implies that the I/Boi/^Edd seems to be 
a better underlying physical parameter than Mbh • 
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